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On the Lagrangian formulation of 
multidimensionally consistent systems 

By Pavlos Xenitidis, Frank Nijhoff & Sarah Lobb 
Department of Applied Mathematics, University of Leeds, Leeds LS2 9JT, UK 

Multidimensional consistency has emerged as a key integrability property for par- 
tial difference equations (PAEs) defined on the "space-time" lattice. It has led, 
among other major insights, to a classification of scalar afhne-linear quadrilateral 
PAEs possessing this property, leading to the so-called ABS list. Recently, a new 
variational principle has been proposed that describes the multidimensional consis- 
tency in terms of discrete Lagrangian multi-forms. This description is based on a 
fundamental and highly nontrivial property of Lagrangians for those integrable lat- 
tice equations, namely the fact that on the solutions of the corresponding PAE the 
Lagrange forms are closed, i.e. they obey a closure relation. Here we extend those 
results to the continuous case: it is known that associated with the integrable PAEs 
there exist systems of PDEs, in fact differential equations with regard to the pa- 
rameters of the lattice as independent variables, which equally possess the property 
of multidimensional consistency. In this paper we establish a universal Lagrange 
structure for afhne-linear quad-lattices alongside a universal Lagrange multi-form 
structure for the corresponding continuous PDEs, and we show that the Lagrange 
forms possess the closure property. 

Keywords: Integrable systems, multidimensional consistency, Lagrangians, 

symmetries 



1. Introduction 

> 

■ The study of integrable systems on the space-time lattice has, over the past decades, 

| developed into a major area of research. The early examples of integrable lattice 

0^ ■ systems, i.e. partial difference equations (PAEs) on two- or higher-dimensional 

grids, go back to the mid 1970s and early 1980s, when the research was focused 
on discretizing known continuous soliton systems (Ablowitz & Ladik 1976, Hirota 
1977, Date et al. 1982,1983, Nijhoff et al. 1983, Quispel et al. 1984). However, 
in recent years the focus has shifted to studying these equations on their own 
merit, since they exhibit the key integrability aspects in a more lucid way than 
their continuous counterparts, and often can be thought of as more generic systems 
than the corresponding differential equations which can be retrieved by appropriate 
continuum limitf[J. 

It has emerged that one of the key integrability aspects of lattice equations 
is the property of multidimensional consistency, i.e., the property that one can 
impose PAEs of the same form (differing only by a choice of the so-called lattice 
parameters) in different pairs of perpendicular directions of a multidimensional 

f It should be pointed out that only very special and non-trivial discretizations of a given 
integrable differential equation will retain the key integrability properties of the latter. 



Article submitted to Royal Society 



TeX Paper 



2 



P. Xenitidis, F. Nijhoff & S. Lobb 



lattice in which the two-dimensional grids are embedded (Nijhoff & Walker 2001, 
Bobenko & Suris 2002). This phenomenon can be understood as the analogue of 
the existence of infinite hierarchies of PDEs associated with an integrable nonlinear 
evolution equation, such as the KdV equation. In some sense multidimensional 
consistency implies that the "integrable system" is not really a single PAE but 
comprises all members of the parameter-family of equations compatible on the 
multidimensional lattice (Nijhoff et al. 2001). Using this property a classification 
of two-dimensional scalar integrable lattice systems was given, in the affine-linear 
case, by Adler, Bobenko and Suris (2003) resulting in what is hereafter referred 
to as the ABS list. Although we shall refer to the equations in the list as the 
"ABS equations", which are given in the Appendix, several examples in the list, 
notably lattice systems of KdV type (associated with discretizations of the KdV and 
related PDEs), were already known in the literature, cf. Nijhoff & Capel (1995) and 
references therein. The most general equation of the ABS list was found by Adler 
(1998) and it contains lattice parameters which are points on an elliptic curve, cf. 
also Adler & Suris (2004). 

Considering the integrable system as a system of many equations imposed si- 
multaneously on one and the same (possibly scalar) dependent variable, raises the 
question whether there exists a canonical description which encodes the whole sys- 
tem of equations. Even though most continuous integrable equations can be cast into 
Lagrangian form (Zakharov & Mikhailov 1980), the situation for discrete equations 
is distinctly more complicated and, in fact, the first Lagrangian for a quadrilat- 
eral equation was given in Capel et al. (1991). Furthermore, the usual least-action 
principle only provides one variational equation per component of the dependent 
variable. Lobb & Nijhoff (2009) provided the first proposal on how to achieve a vari- 
ational description of the multidimensionally consistent system of lattice equations, 
based on the idea of Lagrangian multi-forms. In fact, the proposal is to consider 
action functionals in terms of discrete Lagrange 2-forms Jz? given by 

S = S[u{n);o]=Y J ^= £ -^», (1-1) 

" tJij(n)£tj 

where u(n) is a field configuration of a dependent variable u depending on lattice 
sites labelled by n G T, N of a lattice of arbitrary dimension N > 2 and a denotes 
a quadrilateral surface, composed of a connected configuration of quadrilaterals, 
embedded in the iV-dimensional grid. The least-action principle formulated in Lobb 
& Nijhoff (2009) imposes that S attains a critical point not only with regard to the 
infinitesimal variations of the field variables u(n) for a given chosen surface a , but 
also with regard to the geometry of the independent variables, by imposing surface 
independence of the action for the solutions of the Euler-Lagrange (EL) equations. 
The latter requirement is equivalent to stating that the Lagrangians Jzfy , viewed 
as 2-forms associated with the quadrilaterals that compose the surface, have to 
be closed forms subject to the EL equations. It is this property that makes the 
choice of Lagrangians admissable and which immediately implies multidimensional 
consistency of the system of equations arising from the EL equations on each choice 
of surface a (fixing, if necessary, boundaries). Remarkably, this closure property 
could be established for well-chosen Lagrangians for equations in the ABS list, 
either by direct computation (as was done in Lobb & Nijhoff 2009) or by general 
principles (Bobenko & Suris 2010). Furthermore, Lagrangians possessing the closure 



Article submitted to Royal Society 



On the Lagrangian formulation of multidimensionally consistent systems 3 



property were found in the case of multi-component systems (Lobb & Nijhoff 2010) 
and higher-dimensional systems of KP type (Lobb et al. 2009). 

The phenomenon of multidimensional consistency is not restricted to the dis- 
crete case of lattice equations. In Nijhoff et al. (2000) the novel class of non- 
autonomous PDEs was introduced which are the equations of the same dependent 
variable appearing in the lattice equations but viewed now as a function of a pair of 
(continuous) lattice parameters, denoted by ai,ctj. The so-called generating PDE 
of the KdV hierarchy, i.e. a PDE which generates all the equations in the KdV 
hierarchy by systematic expansion with respect to the lattice parameters, arises as 
the EL equations b u I£ = from a simple Lagrangian given by 

2 d ai ud aj u 2{a i -a j )\ J d a .u o a% u ) 

In Lobb & Nijhoff (2009) it was also shown that this continuous Lagrangian can 
be viewed as a closed 2-form, on solutions of the EL equations, subject to a three- 
dimensional constraint on the solution^. Thus, a similar multidimensional least- 
action principle in terms of a continuous action functional 

S[u(a);a] = fjzf = / da* A day (1.3) 

J" Kj 

which involves continuous Lagrange 2-forms Jz? depending on the fields u(a) and 
on surfaces a with coordinates given by the vector a = (ai, . . . , ckjv), where i,j — 
1,...,N, can be formulated in the continuous case for certain PDEs associated 
with KdV type lattice equations. In order to avoid the difficulty stemming from the 
fact that the scalar Lagrangian (|1.2p . and consequently also the EL equation, is of 
higher order it is natural to formulate the Lagrange structure to that for a system of 
lower-order equations, which actually emerge naturally from the corresponding Lax 
pair. In fact, such systems of PDEs were subsequently established for all equations 
in the ABS list through symmetry analysis (Tsoubelis & Xenitidis 2009) and as a 
specific example of the construction we can mention the following coupled system 
of PDEs, corresponding to HI, H2 and Ql: 



dui Ui - Uj ( , du\ 2 du 

da j ai - otj \ da j J daj 

d 2 u „ Ui — u-j du du ni du n, du „. 

H q — • (1-4) 



daidaj ai — aj dai daj at — <x, da j aj — ai dai 

It was derived in Tsoubelis & Xenitidis (2009) along with its Lax pair, where u, Ui 
and Uj denote a triple of dependent variables. It is worth mentioning that system 
(|1.4[) with (5 = was first given in Nijhoff et al. (2000) in relation to HI and was 
derived from a reduction of the anti-self dual Yang-Mills equations in Tongas et al. 
(2001). In the present paper we will derive the Lagrange structure for the general 
class of systems of the form Jzf (u, Ui, uf, d ai u, d aj u, d ai Uj, d aj Ui) such that the EL 

| In fact, this constraint is equivalent to stating that the solutions obey a KP type of equation 
for the embedding in a higher dimension, which holds true for known infinite families of solutions 
such as soliton solutions. 
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equations, when varying with respect to all three dependent variables, are satisfied 
on solutions of the systems like (11.41) . In this case the Lagrangian reads 

^ du dui du duj fs 2 ( ■) ( Ui ~ u j) 2 \ ® u ® u 

dcti dotj dotj doti \ 1 on ~ ctj J don dotj 

Ui — u-j du iij — Ui du ,„ „. 

- nj J - — + n, — (1.5) 

Cti ~ Ctj OUi dj — di OCtj 

and one can show directly, without assuming further constraints on the higher- 
dimensional embedding, that this Lagrangian obeys the continuous version of the 
closure relation 

d ak J£ij + d ai J£jk + d a .Jz?ki = . (1.6) 

The organization of the paper is as follows. We first give in Section 2 a universal 
description of the Lagrange structure of affine-linear quadrilateral equations. What 
is important here is that the Lagrangians can actually be systematically derived 
from some basic assumptions, such as a 3-point form of the Lagrange function and 
the Kleinian symmetry of the quadrilateral equation. This leads to an integral form 
of the so-called three-leg representation of the quadrilateral equation, cf. Adler et at 
(2003). In Section 3, we specialize to the ABS case, in which we have the symmetries 
of the square, imposing dependence on lattice parameters in a covariant way. We 
prove that the corresponding Lagrangians, in the integral form, obey the closure 
property. In Section 4, we briefly present the construction of the universal system of 
PDEs associated with the ABS list, which was given in Tsoubelis & Xenitidis (2009), 
and we prove the multidimensional consistency of these systems. Finally in Section 
5 we establish the corresponding Lagrangian structures, both for the systems of 
PDEs as well as for the differential-difference equations which are defined in terms 
of the symmetry generators, and show that the universal Lagrange structure for 
the system of PDEs obey the closure relation. Thus, we believe this completes the 
full picture of the Lagrange structure for multidimensionally consistent equations 
associated with the affine-linear scalar quad-lattices. 



2. Lagrangian formulation of 2D quad-lattices 

We first focus on the Lagrangian formulation of affine-linear quadrilateral equations 
on the two-dimensional lattice, and show that under very general assumptions there 
is an almost unique Lagrangian description. The members of the class of equations 
we consider will be denoted by 

Q(u,Ui,Uj,Uij) = 0, (2.1) 

and involve the values of a function u at four neighbouring points of an elementary 
quadrilateral. The labels i, j refer, for the time being, to two fixed directions of a 
1? lattice. 

The general notation we use here, adapted to cover also the situation we will 
encounter in the next section, is as follows. In principle, the unknown function u of 
the equations will, in due course, be considered as functions of an arbitrary number 
of discrete variables, denoted by rij, and there will be an equal number of continuous 
variables on, i.e. u = u(ri\, 11,2, • • • ; Oi\, 012, ■ ■ ■) will be the dependent variable of a 
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system of equations in an arbitrary number of dimensions. In that context, the 
shifted values of u will be denoted by indices, e.g. 

Ui = .Vi(u) := u(...,m + l,...; a), 

u-j = y[ x {u) := u(. ..,rij-l,...;a), 

Uij = .yi-Yjiu) := (...,rii + l,...,nj + l,...;a), 

where S^i is the shift operator in the i direction. Occasionally we will also use the 
difference operator Aj, which is defined by Aj := 5?% ~ 1- The main integrability 
property of the systems considered in the following sections is multidimensional 
consistency. This means that all of them can be extended in a compatible way in 
higher dimensions, according to the following definition: 

Definition 1. Let Q(u,Ui,Uj,Uij) = denote a quadrilateral difference equation, 
to be solved for a dependent variable u. Imposing the dependence of u on a third 
discrete variable nk such that 

Q(u,Uj,u k ,Ujk) = 0, Q(u,u k ,Ui,u ik ) = 0, 

if the three different ways to evaluate the triply shifted value lead to the same 
result, then we call the equation defined by Q multidimensionally consistent. 

However, in the present section we will only make statements about a single 
equation of the form (|2.ip . in terms of two fixed independent variables rii, rij, and 
their corresponding lattice shifts. The following properties of the function Q are 
assumed here: 

1) Q is an affine-linear polynomial depending explicitly on all the arguments, 
that is d u Qd Ui Qd Uj Qd Uij Q £ and 3 2 U Q = d%.Q = d*.Q = dl..Q = 0. 

2) Q is irreducible, meaning that Q cannot be factorized and presented as a 
product of two polynomials; 

3) Q possesses the Kleinian symmetry, i.e., 

Q(u,Ui,Uj,Uij) = eQ(ui,u,Uij,Uj) = a Q(u 3 ,u i3 ,u,Ui) , (2.2) 

where e = ±1 and <r = ±f . 

If Q satisfies Property 1 one can define, in principle, six different bi-quadratic 
polynomials in terms of Q and its derivatives, and four quartic polynomials through 
the application of the following double-sided Wronskian, respectively discriminant 
operator, defined by their actions on an arbitrary polynomial function / = f(x,y) 
in the following ways, eg. Adler et al. (2009): 

®*M ■■= 

Furthermore, the irreducibility of Q guarantees that none of the polynomials 
obtained from Q by acting on it by the operator 2> x ,y for all pairs of arguments 
x, y of Q, is identically zero. The Kleinian symmetry (|2.2I) implies that these bi- 
quadratic polynomials are symmetric and take the same form when we apply the 
operator @ XtV for mutually exclusive pairs of arguments. This is expressed by the 
following Lemma, which can be easily established by exploiting the definitions and 
symmetry relations: 



/ d x f 
d y f d x d y f 



@ x {f) = {d x f) 2 -2fd 2 J. 



(2.3) 
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Lemma 1. Due to the Kleinian symmetry (|2.2[) . the following polynomials assigned 
to an elementary quadrilateral are the same for functions depending on complemen- 
tary pairs of variables: 

1. symmetric bi-quadratic polynomials assigned to the edges of the quadrilateral: 

h(u,Ui) := @ Uj , Uij {Q), h(uj,Uij) := @ u , Ui {Q) , 

h(u,uj) := @ Ut ,u tJ (Q), h(ui,Uij) := @ u , Uj (Q) , (2.4) 

2. symmetric bi-quadratic polynomials assigned to the diagonals of the quadri- 
lateral: 

H( Ul ,u 3 ) := ^u, uy (Q) , H{u,u l3 ) := @ Ui , Uj {Q) ■ (2.5) 

3. quartic polynomials associated with the vertices of the quadrilateral: 

r(u) =: ® Ui (h(u, Ui )) = @ Uj (h(u, Uj )) = $> Uij (H(u, Uij )) . (2.6) 

Furthermore, between the bi-quadratic polynomials the following relations hold on 
solutions of equation \2. 1)) . 

h(u,Ui) h(uj,Uij) = h(u,Uj)h(ui,Uij) = H(u,Uij) H(ui,u 3 ). (2.7) 

With the above definitions we are now in a position to write an integral form of 
the three-leg formula for an afhne-hnear quadrilateral equations of the type (|2.1[) 
subject to the properties l)-3) given above. The usual three-leg formula was given 
in Adler et al. (2003) on a case-by-case basis. In the monograph by Bobenko & 
Suris (2008) a result by Adler is cited as an exercise, of which the following is a 
slight generalization: 

Proposition 2.1. Let Q possess the properties l)-3) given above. Then the equation 

At f Uj dt f u ' j dt 

/ ' (2.8) 



x h(u,t) J y h(u,t) J z ( u ,x,y) H (u,t) 

holds on solutions of Q(u,Ui,Uj,Uij) — 0, where x and y are arbitrary and where 
the function z is defined through the relation Q(u,x,y,z) = 0. 

Proof. The proof is straightforward using implicit differentiation with respect to 
the variables x, y, Ui and Uj. Substituting iiy = z(u,Ui,Uj), since we require rela- 
tion (|2.8[) to hold on solutions of Q(u,Ui,Uj,Uij) = 0, and denoting the resulting 
left-hand side of (I2.8[) by R{u; x,y;Ui,Uj), we can easily prove that R is actually 
independent of x and y, and also of Ui and Uj. To show this it is sufficient to show 
that d x R = through implicit differentiation. Computing the derivative of R with 
respect to x we obtain: 

a« = I ( 2.9) 

dx h{u,x) H(u,z) dx ' 

since z is determined by u, x and y through the equation Q(u, x, y, z) = 0. Implicit 
differentiation of the latter yields 

dz_ _ _ d x Q(u,x,y,z) = _ H(u,z) 
dx d z Q(u,x 7 y, z) h(u,x) ' 
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which holds on solutions of Q(u, x, y, z) — by taking into account the definitions 

h(u,x) = % iZ (Q(u,x,y,z)) , H(u,z) = @ x>y (Q(u, x, y, z)) . 

The combination of (I2.9[) and (|2.10[) leads to the result d x R = 0. In precisely the 
same fashion it can be established that d y R = d Ui R — d Uj R — , where in the 
latter two we make use of the substitution for Uij as mentioned above. Since as 
a consequence R(u;x,y;ui,Uj) does not depend on x,y, we can chose the latter 
arbitrarily, leading to the conclusion that R(u;x,y,Ui,Uj) — R{u;ui,Uj]Ui,Uj) , 
where the latter expression vanishes due to the integration limits and taking into 
account that for x — Ui, y — Uj we have that z = Uij . □ 

Equation f|2 . 8[) suggests that the three leg form can be written as 

F(u, «,-,«,-,««) - F°(u) = 0, (2.11) 

where 

r- dt p at r« dt 

F(u, Ui ,u j>Uij ) := J j^ + J J^ + J Jj^jy ( 2 - 12 ) 

and 

F°(u) := F(u,x,y,z) , (2.13) 

Q{u,x,y,z)— 

where the indefinite integrals denote the primitive (anti-derivative) functions of 
the integrands with respect to the variables Ui, Uj and Uij, respectively, and the 
integration constant in each integral is assumed to be independent of the other 
integration variables. 

Lagrangians for the generalized three-leg formula (|2.8p can be obtained by simple 
integration over the remaining variable u, but there is a stronger result that holds. 
In fact, it turns out that under the assumption of the three-point form of the 
Lagrange function, as was done in Lobb & Nijhoff (2009), one can actually derive 
in an almost unique way the form of the Lagrange function for the affme-linear 
quadrilateral equations. This is expressed in 

Theorem 2.2. Let Q(u, Ui, Uj, Uij) be an irreducible, affine linear and Kleinian- 
symmetric polynomial. Then, up to a freedom of having a multiplicative constant, 
and an additive constant, a three-point Lagrangian of the form J£(u,Ui, Uj) of the 
equation Q = 0, is given by: 



^£{u, Ui, Uj 



dsdt m dsdt r-n dsdt 



h{s,t) J J h(s,t) J J H(s,t) 
F°(s)ds. (2.14) 



Integrals denote anti- derivatives of the integrands which are determined up to con- 
stants, i.e. independent of u, Ui andiLj. 

Proof. Let us suppose that we are looking for a function Ui, Uj) such that 

<5j£f d / \ 

~~6u ^ Q^\^( u ' u i' u j) + ^?( u -ii u , u -i,j) + ^'(u-j,Ui^ j ,u)J = (2.15) 
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on solutions of the equation Q = 0. 

First we differentiate the above equation with respect to it_i by taking into 
account that U-ij is determined by the equation Q(u—i,u,U—i t j,Uj) = 0. This 
yields 

Qn - ■ 

{u~i,u, u-ij) + du^ , d u J£(u-i, u, U-ij) _ ~*' 3 = 0, (2.16) 

OU-i 

where 

du-i,j _ d u _ i Q(u- i ,u,u- i j,Uj) _ H(u,u-ij) 
d u _ i:l Q(u-i,u,U-i,j,Uj) h(u,u-i) 

on solutions of equation Q(u—i,u,U—ij,Uj) = 0. Thus, substituting the latter into 
(j2.16[) and shifting forward in the i direction the resulting relation, one arrives at 

h(u,Ui) d u d Ui J£(u,Ui,Uj) = H(ui,Uj) d Ui d Uj ^f(u,Ui,Uj) . (2.17) 

In the same fashion, the derivative of (|2.15[) with respect to U-j leads to 

h(u,Uj) d u d Uj 3f(u,Ui,Uj) = H(ui,Uj)d Ui d Ui J£{u,Ui,Uj) . (2.18) 

A third equation follows from the compatibility of (|2. 1712.181) . namely 

h(u,Ui)d u d Ui 3?(u,Ui,Uj) = h(u,uj) d u d Uj ^f(u,Ui,Uj) . (2.19) 

The solution of the system of equations (|2.17fj2~T9|) is 

/ rr* dsdt rr^dsdt r Ui f u * dsdt \ 

(2.20) 

where c is a constant, $, ^ and X are arbitrary functions of their arguments and 
double integrals denote anti-derivatives. 

To determine the unspecified functions involved in Jzf, we write out explicitly 
the Euler-Lagrange equations (|2.15p using (|2.20p and taking into account (|2.11l 
|2~T3|) . This leads to 2cF°(u) + + *'(it) + X'(tt) = 0, where prime denotes 

differentiation with respect to u. From the latter equation, one of the unspecified 
functions can be defined in terms of the remaining ones. We choose to solve it for 
$(w) leading to 

f u 

$(u) = - X(tt) -2c F°(s)ds, 

which is defined up to an additive constant. 

Substituting back into (|2.20l) . it is not difficult to see that the terms involving 
functions ^ and X can be written as a total difference, and since their variational 
derivative does not contribute to the EL equations, they can be omitted from Jz? . 
Hence, up to a multiplicative constant and an additive constant, Jzf is given by 
(EH- □ 

In all of the above formulas, primitives are determined up to constants. One can 
fix this freedom during the derivation of the Lagrangian by choosing an arbitrary 
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function u° along with its shifts as lower limits in the integrals involved in (|2.20p . 
The derivation determines the remaining terms in the latter leading to 



_£f (u, Ui, uj) 



dsdt r r j dsdt r ui r* dsdt 



h(s,t) J u oJ u o h(s,t) Ju Ju° H(s,t) 



H(s,t) J u o J u o H(s,i) 



,(2.21) 



where w and z are solutions of the equations Q(u° , s, w, ity) = and Q(u , z, t, = 
0, respectively. 

Lagrangian (|2. 14[) possesses the Klein symmetry in the sense that three different 
Lagrangians can be derived from this one by applying the interchanges implied by 
the Klein symmetry. 

To make this statement clear, we apply the changes (|2.2[) to (|2.14l) . The first 
interchange, i.e. (u, Ui, Uj, utj) — > (u%, u, utj, Uj), which leaves invariant the equation 
Q = 0, when applied to «5f leads to «£?(uj, which is another Lagrangian of 

equation Q = 0. Additionally, applying the changes (u, Ui, Uj, Uij) —> (uj,Uij, itj, it) 
to (|2.14p one arrives at another Lagrangian, namely J£(uj,Uij,u). Finally, combin- 
ing these two symmetries one can derive a third Lagrangian which has the form 

Uj,Ui). 

In Bobenko & Suris (2010) the proof of the closure relation makes use of a fun- 
damental property of the Lagrangian, which in the context of the present approach 
can be generalized and understood as the following 

Proposition 2.3. All the above Langrangians are defined up to an additive constant 
and, by choosing it appropriately, their sum is zero on solutions of the equation 
Q = 0, i.e. 

J£(u,Ui,Uj) + Jf{ui,u,Uij) + J£{uij,Uj,Ui) + Jf(uj,Uij,u) = 0. (2.22) 

Proof. To prove that the sum is a constant on solutions of Q = 0, it is sufficient to 
prove that its derivatives with respect to u, Ui and Uj vanish. Indeed, differentiating 
it with respect to u or the single shifted values of u and taking into account that 
depends on u, Ui and Uj through the equation Q = 0, one easily verifies that 
the sum of the four Lagrangians is constant on solutions of Q = 0. □ 



3. Universal Lagrangian formulation of the ABS equations 

and Closure 

Among the members of the class of equations studied in the previous section are 
the ones possessing the symmetries of the square, or D4 symmetry. These equations 
may depend on two lattice parameters, denoted by at and a*j respectively, and their 
defining polynomial will be denoted by Qij(u, Ui, Uj,mj), or simply by Qij, in order 
to make this dependence explicit. 

Since the lattice parameters are assigned to the edges of a plaquette, the sym- 
metries of the square imply that the interchange of Ui and Uj must be followed by 
a mutual interchange of the lattice parameters, i.e. 

Qij(u,Ui,uj,Uij) = eQji(u,Uj,Ui,Uij) . (3.1) 
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The ABS equations, which are given in | Appendix A[ belong to this subclass and 
depend explicitly on the lattice parameters. They are multidimensionally consistent, 
cf. Definition [1] and this property depends essentially on the appearance of the 
lattice parameters in them. These two characteristics of the ABS equations have 
the following consequences. First of all, the polynomials h assigned to the edges can 
be factorized as 

h(u,Ui) = Kijhi(u,Ui) , h(uj,Uij) = Kijhi(uj,Uij), 

h(u,Uj) = Kjihj(u,Uj), h{ui,Uij) = Kjihj(ui,Uij). (3.2) 

The function Kij — K(a,i,a.j) is antisymmetric, i.e. Kji = —Kij, and hi(x,y) is 
a symmetric bi-quadratic polynomial of x and y and the index i denotes that it 
depends only on the lattice parameter oti. 

Moreover, polynomial H(x,y) is symmetric and bi-quadratic with respect to x 
and y, and symmetric with respect to the lattice parameters. It is more convenient 
for the analysis of the ABS equations in the following sections to introduce the 
function 

Hyfry) := (3.3) 

which is antisymmetric with respect to the interchange of indices : Hij = —Hji. 
Again, the indices denote the dependence on the lattice parameters and not shifts 
in the corresponding lattice directions. Functions hi, H and Kij can be found also 
in | Appendix A| 

Finally, the factorization of polynomials h implies that the three leg form for 
these equations can be written as 

: dt _ n dt + r« tit _ Q 
hi{u,t) J y hj(u,t) J z{ 

u.x.y 

where x, y are arbitrary and z is any solution of Qij(u, x,y, z) = 0. 

In this case, the function F (]2 . 12[) is antisymmetric with respect to the inter- 
changes of the indices and is given by 

f u ' dt f Uj dt f Uii dt 

Fij ee F(u,Ui,Uj,Uij;a>i,aj) := J ^ - J + J • (3-5) 

Similarly, the function F° (|2 . 1 3[) takes the following form 

F°Au) := F(u,x,y,z;ai,aj) . (3.6) 

Qij (u,x,y,z)=0 

Having introduced through Proposition l2.2l a Lagrangian for the ABS equations, 
we now proceed and prove that it obeys the closure property. 

Theorem 3.1. The following function 

dsdt ff"' dsdt r Ui r Uj dsdt 



hi(s,t) J J hj(s,t) J J H rj (s,t) 
F?At)dt- / FfJt)dt (3.7) 
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defines a Lagrangian for the ABS eguation Qij = 0, meaning that the EL equations 
are satisfied whenever Qij(u,v,i,uj,Uij) — 0. Moreover, the Lagrangians J£jt 
and ££n satisfy the closure relation 

A^jt + A^ u + AtJfij = , (3.8) 

which holds on solutions of the ABS equation Qij = 0, Qjt — and Qu = 0, 
and hence, in view of skew symmetry of the Lagrangian functions Jz?ij(u, Ui, Uj) — 
i), it defines a closed discrete 2- form. 

For the proof of the closure relation, we will need the following 

Lemma 3.2. Let hi(u,x) be any of the polynomials related to the ABS equation 
Qij = 0. Then, the relation 

2d CH h i (u,x) + 9 UiX {hi{u,x)) = (3.9) 

holds identically. 

Similarly, if Hij(u, x) is any of the diagonal polynomials corresponding to the 
ABS equation Qij — 0, then the following identities hold 

2d ai H ij (u,x)-@ u , a: (H ij (u,x)) = 0, 2d aj H ij {u,x)+9 UtX {H ij {u,x)) = 0. (3.10) 

The identity for hi follows from the symmetry analysis of the ABS equations, 
(Xenitidis 2009, Mikhailov et al. 2010), and its derivation is given in |Appendix B 
Actually, this relation expresses the fact that the lattice parameters play the role of 
the continuous "master symmetry" parameters. The same relation can be found in 
the context of integrable chains of differential-difference equations and their master 
symmetries (Yamilov 2006). 

The relation for ify can be proven either directly, or, in the case of Q1-Q4, 
using the observation that Hij(u, x) = —h(u, x, a; L — ay). 

We now proceed with the proof of Theorem 13. II 

Proof. Since the ABS equations belong to a subclass of the family of equations 
studied in the previous section, the Lagrangian (|3.7p easily follows from (|2.20l) by 
assuming it to be symmetric with respect to the interchange of indices, choosing 
the multiplicative constant to be Kij and using relations (|3.2l 13. 3p . 

Let us now denote the left hand side of (|3.8p by ^ . To prove that the closure 
relation holds on solutions of equations Qij — 0, Qje = and Qu — 0, it is sufficient 
to prove that ^ is independent of u and its single shifts, as well as the lattice 
parameters. Observing that the closure relation is symmetric with respect to any 
permutation of the indices, it is sufficient to show that d u ^ — — d ai c € = 0. 

Actually, it can be proven that all the first order derivatives of ^ vanish on 
solutions of Qij = 0, Qji = and Qu — 0. First of all, writing ^ out explicitly, it 
follows that it is independent of u. Now, on the one hand, the derivative of ^ with 
respect to Ui reads as follows 

d u ^ = -T t] {u % )-T H {ui), (3.11) 

where 

Tij(ui) := Fij(u,Ui,Uj,Uij) - F°j(in). 
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Obviously, in view of Qij = and Qa — 0, relation (|3.11l) implies that d u ^ = 0. 
On the other hand, the derivative of ^€ with respect to Uy has the form 

d Uij V = y^T^Uj)) + ^(T«(ui)) = 0, (3.12) 

where the last equality holds on Qji = and Qu = and their shifts. Thus, 
due to its symmetry, is independent of u and all of its shifts. 

The next step is to prove that the closure relation is independent of the lattice 
parameters. Considering the derivative of ^ with respect to cti and using Lemma 
(13.21) . one arrives at 

d <£ — - in f Hu('u,v,i j e)Hei(ii,jt,Uij)\ dl 

where 



2 \Hij(u,Uij)Hij(uij,Uje) J da 



(pUi rUj \ / put rUi 

J Ffj(8)d8 + j F^dsj+AJJ F«(,s)d s + y F^(s)ds 

It is not difficult to verify case by case that 
Hu{u,uu)Hu{ujz,Uij) - Hij(u,Uij)Hij{uu,Uji) 



}ij=o,Qji=o,Q u =n 

which is equivalent to the affine linear equation P(u,Uij,Uji,uu) = relating the 
four points and Uit. 

Finally, the definition (|3 . 6[) of i 7 * and relations!}] 

duj-j _ -Hjjfaujj) j / hj{u,Uj) \ dujj _ -Hj^u^Uj) ( hj(uj,Uij) 
dcti 2 " \Hij(u,Uij) J ' don 2 " J \Hij(ui,Uj) 

(3.13) 

imply that d a J = 0. Combining all these it follows that = 0. 

Hence, the closure relation is constant on solutions of the corresponding ABS 
equation. We will show now that this constant is zero. The left hand side of the 
closure relation, when is evaluated on solutions of the corresponding ABS equation, 
becomes 



Hij(s,t) J J Hij(s,t) J J hi(s,i) J J hi(s,t) 

+ terms obtained by cycling permuting indices . (3-14) 

This formula suggests the existence of the solution^ 

Uij = ut, Uji = Ui, ua = Uj , (3.15) 

which, when combined with relation (|3.14l) . leads to = 0, implying that van- 
ishes on any solution of the corresponding ABS equation. 

f These relations hold on solutions of the equation Qij = and can be verified case by case 
since their proof is lengthy and is omitted here. 

| This solution was considered also in Bobcnko & Suris (2010). 
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The particular solution (13. 15)) is equivalent to the system 

Qij(u,Ui,Uj,ui) = 0, Qjt(u,Uj,ue,Ui) = 0, Qii(u,ui,Ui,Uj) = 0. 

To prove that such a solution exists, it is sufficient to prove that two of the equa- 
tions imply the third one. Indeed, writing the two first equations in terms of the 
polynomials h, i.e. 

hi{u,Ui)hi{uj,ui) = hj(u,Uj)hj(ui,Ui), hj(u,Uj)hj(ue,Ui) = ht(u,ui)ht{uj,Ui), 

the elimination of polynomials hj between them yields 

hi(u,Ui)hi(uj,ue) = h e (u,ue) h e (uj,Ui) , 

which is equivalent to the third equation of the above system. □ 

Remark 3.1. There exist deformations of the H equations in the ABS list, which 
were first presented in Adlcr et al. (2009), and further studied in Xenitidis & Papa- 
georgiou (2009). The defining functions of these equations possess the basic proper- 
ties of affine linearity and irreducibility, and the symmetries of the rhombus. They 
are defined on a black-white lattice which implies that they can be written in non- 
autonomous form. In this formulation, the polynomials assigned to the edges have 
the same form and can be factorized like the polynomials (|3.2p corresponding to 
the ABS equations. On the other hand, there are two diagonal polynomials Hi and 
H2 related to each other by shifts, e.g. S^iHi(ui,Uj) — #2 («« , ) , cf. Xenitidis 
& Papageorgiou (2009). Actually, from Theorem 13.11 follows a Lagrangian for the 
deformed H equations which satisfies the closure relation as well. One has to use 
the polynomials h corresponding to the deformed equations and replace Hij(u, uij) 
and Hij{ui,Uj) by HnjiuijUj) and H2ij{u,Uij), respectively, in Q3.7p . □ 

Remark 3.2. The lift of the ABS equation Q,j = to a two-field discrete system 
was suggested by Papageorgiou & Tongas (2009) in the construction of Yang-Baxter 
maps. Denoting by u and v the fields involved in this system, the latter has the 
following form 

Qijlu^v^u^u^) = 0, Qij(v,Vi,Uj,Vij) = 0. (3.16) 

A Lagrangian for the above system follows actually from p.7p and it has the fol- 
lowing form 



JZ?(M, V,Vi,Uj) 



dsdt m dsdt dsdt 



hi(s,t) J J hj(s,t) J J Hij(s,t) 
h F°(s)ds - fk(s)d S , 



where F?- is defined in (|3.6p . It is a straightforward calculation to prove that the 
variational derivatives of Jz? (u, v, Vi, Uj) with respect to u and v vanish on solutions 
of Qij(v-j ,Vi.-j,u,Vi) = and Qij(v,-i,V,U-i t j,Uj) — 0, respectively. □ 

f Here we have used the fact that polynomials hi are symmetric, i.e. hi(x,y) = hi(y,x). 
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4. Multidimensionally consistent continuous systems from 

symmetry reductions 

All of the ABS equations admit a pair of extended generalized symmetries which 
play the role of master symmetries for the corresponding generalized symmetries 
(Rasin & Hydon 2008, Xenitidis 2009). This pair of extended symmetries can be 
given as the differential operators 

Vi = m \ % U,Ul ' - \d Ui hi{u,m) J d u - d ai . (4.1) 
\Ui — U-i 2 J 

Since the equations under consideration depend explicitly on the lattice param- 
eters, the same holds for their solutions. Thus, it is interesting to consider solutions 
of the ABS equations remaining invariant under the action of the extended sym- 
metries. Such solutions satisfy the discrete equation under consideration and the 
system of differential-difference equations 

du ( hi(u,Ui) 1 \ 

-d Ui hi{u,Ui) =0, (4.2) 



From this overdetermined system Tsoubelis & Xenitidis (2009) derived the following 
system of partial differential equations 

dui du _ „ , , „. 

dai hi da. 



x 3 J WLX J 

duj Hji du 



rii ViHj, , (4.4) 



dai hi dai 

d 2 u du du „ du _ <9it _, , , 

ir~ + B a s~ + + r y > ( 4 - 5 ) 



daidaj ctai Sa., <9aj 

where the arguments of hi(u,Ui), hj(u,Uj) and Hij{ui,Uj) have been omitted and 
n,i, rij were scaled by 2. Moreover, 

V t A := ^ - lnfc) X = ^- f J V (4.6) 

and the coefficients in the last equation of the system are given by the following 
relations. 

Aij := d u In hi - j-^i H ij = d u In hj + ^- Vj ffy , (4.7) 

/<',, : -nj V.Vj //,, = rij (-hj d u Vjl + ^-V ,;// , Vj //,,.;) , (4.8) 
r^j := -n l n j V t (hjS7 i 'S7 j Hi-j) = mrijV 'j(hjV jV 'j ify) , (4.9) 
In what follows, system (|4. 31 14.9ft will be denoted as Sy. 

Remark 4.1. The derivation of Sy implies that u; and Uj are related to u by 
shifts. Actually, this relation can be forgotten and all the functions involved in Sy 
can be treated as independent, which is the point of view we adopt from now on: we 
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consider Ey as a closed-form coupled system for the three functions u, Ui and Uj. □ 

Remark 4.2. The two equivalent forms at the definitions (|4.7p - (|4.9l) of the coeffi- 
cients Aij , Bij and I\j , respectively, follow from the two different ways which can 
be used to derive the second equation of the system and the fact that the master 
symmetries commute (Xenitidis 2009). The latter implies that these equalities are 
identities. In fact, the basic identity is relation (j4.7[) from which (|4.8l4.9p follow. 
Specifically, the V^-derivative of (14.71) leads to (|4.8p . On the other hand, multiplying 
the latter by hi and taking the V^-derivative of the resulting equation, one arrives 
at the third identity (|4~9l) . □ 

Some useful general properties of the operator Vi defined in (I4.6P can be derived 
directly from its definition and are given in the following statement. 

Proposition 4.1. Let hi(u,Ui), hj(u,Uj) be any pair of smooth functions and 
be defined as in |^.6p . Then the following identities hold. 

i) Vi/ii = identically; 

ii) [Vi , Vj] X — 0, for any function X ; 

Hi) d u ViVjY + (d u d Ui hihi)VjY + (d u d Uj lnhj)ViY = 0, for any function Y in- 
dependent of u. 

An important property of E,-j is its multidimensional consistency, i.e. the com- 
patibility conditions among the equations of systems Ey, Hji and E# lead to no 
further restrictions. To make more precise what we mean by a multidimensionally 
consistent system of PDEs, we give the following definition in a way parallel to 
Definition Q] of Section [5] for the discrete case. 

Definition 2. Let F(u,d ai u,d aj u,d ai d aj u) = denote a system of differential 
equations. Imposing the dependence of u on a third continuous variable ak such 
that 

F(u,d aj u,d ak u,d a] d ak u) = 0, F(u, d ak u, <9 Qi u, d ak d at u) = 0, 

if the three different ways to evaluate d ai d aj d ak u lead to the same result, then we 
call the system of PDEs defined by F multidimensionally consistent. 

The multidimensional consistency of E^ can be checked by straightforward cal- 
culations and is equivalent to the identity 



which can be checked case by case for all the ABS equations and their corresponding 
functions ffy. 

To see how (|4.10p expresses the multidimensional consistency of the ABS equa- 
tions and their continuous counterparts Ey, consider the compatibility condition 



which must hold on the solutions of Ey, T,jg and E^. The latter relation, when it 
is written out explicitly using the equations constituting these systems, leads to 



(Hj£d ue Hgi — Hgid ue Hje) + {Hud Ui Hij — Hijd Ui Hu) 
+ (Hijd Uj Hjg — HjidujHij) = , 



(4.10) 





4 



hjhiVjVe = , 
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which holds in view of (|4. 10[) . In the same fashion, the three different ways to evalu- 
ate d ai d aj d ak u and their consistency lead to similar expressions which are satisfied 
whenever = 0. In the following section where Lagrangian functions will be 
constructed for the system E^- we will see that the latter identity also ensures that 
these Lagrangians satisfy the closure relation. 

Remark 4.3. As for the discrete ABS equations, where a Lax pair follows from 
multidimensional consistency (Bobenko & Suris 2002, Nijhoff 2002), the multidi- 
mensional consistency of E^ can be used to derive an auto-Backlund transformation 
and a Lax pair for this system (Tsoubelis & Xenitidis 2009). To present them, it is 
more convenient to introduce the following notation: 

~ H(ui, u, ai, A) 

-Hi : = ; rr , 

K(a i7 A) 

i.e. Hij with Uj and aj replaced by u and A, respectively. Now, the auto-Backlund 
transformation can be written as 



du Hi du „ - du Hj d 



XL - 

o — = ~r- « l~ n i Vi-ffi , -7: — = -jf- 1- n j VjHj , 

oai hi oai daj hj oaj 

Q(u,Ui,u,v,i; ai, A) = 0, Q(u,Uj,u,Uj;a.j,X) = 0, 

where the last two equations coincide with the discrete equation related to Ejj. 
Moreover, a Lax pair for Ey has the following form 

— = Mi $ , — = Mj $ 
eta, Oai 



where 



-I A' — I A" 
4 I A' 



with 



. Hi du - 

Ai := — 1- rii ViHi 

hi da t 



and prime denoting differentiation with respect to u. Also, in the matrices Mj, the 
function Ai and its derivatives are evaluated at u — 0. □ 



5. Lagrangian formulation of the differential-difference and 

PDE system 

In this section, we will show that there exist Lagrange structures for the differential- 
difference equations (|4.2j) and for the PDE system E^ , where in addition the latter 
satisfies a closure relation. We will start by presenting the Lagrangian for the system 
defined by the differential-difference equation 



Proposition 1. The following action functional 



dx 



o hi(u-i,x) 



ln(iii - u-i) - - \n(hi(u, m)) 



du'l 
dai 



da: 



~ln(hi(u,uD) 



u o hi(x 7 u°) 
dai 



(5.1) 
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forms an action for the differential- difference system (I4.2[) in the sense that the EL 
equations are satisfied if u — u(ni,cn) is a solution of the latter. 



Proof. The proof is by direct computation, and uses crucially the relation (|3.9[) , as 
well as the identity 

hi(u,Ui) 1 hi(u,u-i) 1„ , . 

-^Ouihi(u,Ui) = h -Ou^hi^u-i) , 

Ui — u_i 2 m — I 

which holds for general symmetric bi-quadratics hi. Thus, we obtain the EL equa- 
tions in the following form: 



dSe 1 du- 4 1 8 



du hi(u,u-i) dcii hi(u,Ui) dcti 

. rii + 1 , rii - 1 1 rii + 1 . 1 Hj-1 s 

-d u hi{u,Ui) - - — -d u hi(u,u_i) 



u-Uij u-u-i^i 2hi(u,Ui) ' 2hi(u,u-i 

where Jzf denotes the expression between curly brackets on the right-hand side of 
(15. ip , which is a combination of two copies of the differential-difference equation 
(14. 2[) one shifted forward, and one shifted backward in the lattice variable rii. d 

The multidimensional consistency of can be established on the level of its 
Lagrangian formulation by proving the closure relation for the corresponding La- 
grangian 2-form structure in accordance with its discrete counterpart. 

Theorem 5.1. A Lagrangian of system Ejj is given by 

= i_d^(M _ i_ dv_ duj_ _ jL^d^dv_ _ v v . H . (5 2) 

VJ hi dcti dctj hj dctj dcti hi hj dcti dotj 1 3 1 3 13 

1 „ „ „ , , \ 3ii / 1 „ „ „ , , \ 9u 

H 



— Vj Hij + d u In hj j Q— + n» I — V; Hji + d u In hi J 



and is antisymmetric: S£ji = — J£ij . 

Moreover, Lagrangians S£ij , J£ji and S£u satisfy the closure relation 

T> ai ^ji + T> aj 2ti + Do,^- = (5.3) 
on solutions of systems Sjj, Tiji and E^. 

Proof. The variational derivatives of Jzfy with respect to Ui and Uj yield the sec- 
ond equation of Ey, while its variational derivative with respect to u leads to a 
linear combination of the two first equations of Ey . Indeed, consider the variational 
derivative of Jzfjj with respect to i.e. 

SJCij _ d ( 1 du \ d&ij 



Sui dctj \ hi dcti J d u i 

Writing out explicitly the above relation, the terms involving the derivative of u, 
with respect to a.j cancel out. Factorizing the remaining terms, one arrives at 

S^ij 1 / d 2 u du du du du \ 

Sm hi \dctidctj 13 dcti dctj K 'dcti 31 'dctj 13 ) 
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which is identically zero in view of the third equation of . 

Similarly, the variational derivative of the Lagrangian with respect to Uj takes 
the form 

53%j _ -1 / d 2 u _ A _du_du_ _ B 9u du _ , _ () 

5uj hj \daidctj 13 dai daj 13 dai 3 1 daj 13 

For the variational derivative with respect to u, we have 



SJifij d ( 1 / dui du 

5u dai \ hi \ dai hi da "'"/"<.;( .,<„.<.,,, 



— rij V 7 - Hij ) — rij d u In hj 

3 '"3 w ^3 

+ W~ \t~ (w^~ + IT W~ ~ n i^i H ij^) + m du lnh\ - ^ . 

daj I hj \dai hi dai J J du 

The terms in the parentheses are the two first equations of the system, thus the 
above equation becomes 

<L§fj d d dJ£i ■ 

= -— (-rijd u lnhj) + - — (n;9„ln^) — ^ . 

Su dai daj du 

Writing out explicitly the right hand side of the last relation, we use the two first 
equations of the system to substitute the derivatives of Ui and Uj. After these 
substitutions, the terms involving the derivatives of u cancel out and we arrive at 



i fSuVjVj Hij + (d u d Ui Inftj)Vj H l3 + (d u d Uj \nhj)Vi Hi 



S^ij 

Su 

Using the third identity of Proposition l4.il one easily concludes that <5 u Jz?ij = 0. 
Denoting by Si the left hand side of f|5 . 3[) . i.e. 



we will prove that S% = on solutions of Ey, Hji and 

First, we work out explicitly ffl to show that the second order derivatives of Uj, 
Uj and ui cancel out. Then, we factorize the remaining terms with respect to the 
derivatives of it, which leads to the following expression 

d 2 u ( 1 duj 1 dui 4 du A du „ 



Aw- h A, 



where 



oaidaj \hj oae hi oae dai dai 

+ similar terms obtained by cycling permuting indices (i, j, £) 
+ terms involving only first order derivatives of u , (5-4) 



H ( 1 



The coefficient of d ai d aj u in (|5.4I) vanishes by taking into account the first 
order equations of T,ji and £«. Due to the symmetry of the coefficients of the 
remaining second order derivatives of u also vanish. 
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On the other hand, the terms in (15.41) involving only first order derivatives of u 
take the following form: 

1 du du du 



hihjhi dcti dctj dat 

no du du rii du du m du du , 

- V* + ~ 7 — — Vi + rr4- — — Vj T ijt 



hihj dcti dctj hjhe dcxj dai hihg dai dcti 

ujm du mm du mnj du \ 

—r — q — VjV£ H — - — V^Vi H — - — ViVj ±ije. 

hi dai /ij ooj he. don J 

The multidimensional consistency of the continuous systems, i.e. identity (|4.10l) . 
implies that the above expression vanishes on the solutions of systems , Hj£ and 
T,u and, subsequently, the closure relation (15.31) is satisfied. □ 

P.X. is supported by the Newton International Fellowship grant NF082473 entitled "Sym- 
metries and integrability of lattice equations and related partial differential equations", 
which is run by The British Academy, The Royal Academy of Engineering and The Royal 
Society. S.L. was supported by the UK Engineering and Physical Sciences Research Coun- 
cil (EPSRC). 

Appendix A. The list of the ABS equations and their 

characteristic polynomials 

The following constitutes the list of scalar affine-linear multidimensionally consis- 
tent quadrilateral lattice equations up to Mobius transformations that was estab- 
lished by Adler, Bobenko & Suris (2003). The equations of the form Al, A2 of 
their paper are omitted as they are related to other members of the list by point 
transformations. The form of Q4 given below was established by Hietarinta (2005), 
in which sn denotes the Jacobi elliptic function sn(x\k) with modulus k. 

(u — Uij) {ui — Uj) — on + ctj = (HI) 

(u — Uij)(ui — Uj) + {ctj — oii)[u + Ui + Uj + u^) — af + oq = (H2) 

e -a,/2( UUi + u . u ..j _ e -«i/2( Wj + u . u ..} + 8( e -Ci _ e-oy) = Q (H3) 

Oi{u — Uj)(ui — Uij) — ctj{u — Ui)(uj — Uij) + 5 2 aiaj(ai — otj) = (Ql) 

Cti{u - Uj)(Ui - Uij) - Ctj{u - Ui)(uj - Uij) + 

otiOLjiui — otj)(u + Ui + Uj + Uij) — aiCtjipti — Cij)(cq — otiOLj + aj) — (Q2) 
(e - " 3 — e~ ai )(uuij + UiUj) + e~ a ^ 2 (e~ ai — l)(mtj + UjUij) 

_ e - ai{e - aj _ 1)(wi . + u . Uij) _ s^^-^n^-r,^-!) = (Q3) 

sn(ai)(wwi + UjUij) — sn(aj)(uuj + UiUij) — sn(a; — ctj)(uuij + UiUj)+ 
ksn(cti) sn(a.,) sn(aj — aj) (1 + uUiUjUij) = (Q4) 
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HI : 


hi = 1 




h 'ij 


— a,j — ai 


H2 : 


hi = 2(u- 


-tii + a t ) 


Kjj 


— aj — ai 


H3 : 


hi = uui - 


be -o<i/2 § 


Kij 





The formulae for the canonical bi-quadratics h and H , and the corresponding 
discriminant curves, for the ABS equations are collected in the following lists. 

H = (u,i -Uj) 2 
H = (it; - u 3 ) 2 - (ai - aj) 2 
H = e-^+ a ^/ 2 (u 2 +u 2 ) 
— (e~ Ql + e~ a] )u t Uj 

Ql: hi = ((u - Ui) 2 - a 2 S 2 )/ai , = aiaj(aj - a^ 

H = a t aj ((ui - Uj) 2 - (a, - otj) 2 8 2 ) 
Q2 : hi — ((u - Ui) 2 -2af(u + Ui) + a\)ja l , K,j — a t a 3 (ay - ai) 

H = a t a 3 ((14 - Uj) 2 - 2(aj - ay) 2 (ui + Uj) + (a l - ay) 4 ) 

QO. u _ cxpt-ai/2) / 2 | 2\ l+cxp(-a t ) r 2 l-cxp(- Qi ) 
■ — l-exp(-a<) \ U ^ U i> l-cxp(- Ql ) 1 4cxp(-a;/2) 

= [exp(— a») — 1] [exp(-ay) - 1] [exp(-aj) - exp(-ay)] 

H = (c- a - - l)(e- a i - l)(e-( ai+a ^/ 2 (u 2 + u 2 ) - (c~ Q * +e~^)u lUj 

+ I S 2 e ( Q '+^)/ 2 (e" Q * - e-°i ) 2 ) 

Q4 : hi = -fcsu(af)(l + u 2 uf) + (u 2 + u 2 - 2cn(a!t)dn(a:j)mij) 

= sn(aj)sn(ay)sn(ay — ai) 

H = sn(aj)sn(aj) {-fcsn 2 (ai - aj)(l + u 2 u 2 ) + u 2 + u 2 } + 

(sn 2 (ai — ay) — sn 2 («i) — sn 2 (ay) + fc 2 sn 2 (ai — ay)sn 2 (ai)sn 2 (ay)) UiUj 





HI 


H2 


H3 Ql 


Q2 Q3 Q4 


r(u) 





4 


u 2 4S 2 


16m u 2 - 6 2 4 (ku 4 - (1 + fc 2 )u 2 + k) 



Appendix B. Proof of Lemma 13.21 

As was proven in Tongas et al. (2007), all of the ABS equations admit three point 
generalized symmetries which have the following form 

= P(hi(u,Ui)) d u := j ( — - \ du ^j d u ■ 

Moreover, it was shown that any other three-point generalized symmetry K is 
necessarily of the form K — ^a(n) + cj)(n, m,it)| d u , where the functions 
a(n), (j)(n,m,u) are determined by the solutions of a linear equation. Additionally, 
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it was proven that all of the ABS equations admit extended generalized symmetries, 
namely Vi given in (|4.1[) . 

On the other hand, it was shown in Mikhailov et al. (2010) that a higher order 
generalized symmetry is given by 



K 



(2) 



hj(u, Ui)hj(u, u_j) 



i 



9u 



The commutator of the extended symmetry V. with the first generalized sym- 



metry yields 



ViyK™] = -K\*> + P ( d a Mu,Ui) + - 



-(2) 



1 



which is another symmetry of the equation under consideration. Since is a 

symmetry of the equation, the same must hold for the remaining part of the above 
commutator. This means that the latter has to be of the general form K which 
implies 



d ai hi(u, + 



hi d u hi 
d Ui hi d u d Ui hi 



/z(ai) hi(u, u^ . 



But, for all of the ABS equations, it turns out that /i = 0. 
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